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BLOCK PROGRAMMING 


A Valuable Aid to Simulation of Lumped-Parameter Physical Systems 


The analog computer, or differential analyzer, isa 
powerful tool for obtaining dynamic solutions to differ- 
ential equations. Generally the equations to be sol- 
ved relate to a physical system of particular interest 
to the computer programmer, in which case the pro- 
grammer wishes to simulate his physical system by 
means of an electrical network whose defining equa- 
tions are analogous (hence the name analog computer) 
to the subject physical system. It is desirable then, 
to arrange the computer mechanizationto have a one- 
to-one correspondence between the problem board and 
the physical system, so that the programmer ‘‘sees" 
his system rather thanan abstract electrical network. 


The following discussion will illustrate a method of 
programming in which the physical elements are 
treated as blocks rather than as sets of differential 
equations. Certain useful basic building blocks will 
be developed, and methods for combining the basic 
blocks to form any desired physical system (simula- 
tion) will be discussed. Simulations generated by the 
methods to be described will have the advantage that 
any or all physical parameters appear as single con- 
trols in the computer so that any parameter may be 
changed without affecting any other setting. In addi- 
tion, any or all variables appear as computer outputs 
and may be individually monitored. 


A lumped-parameter system, by definition, consists 
of a number of interconnected discrete elements whose 
individual or collective responses to impressed forces 
or stimuli are of analytical interest. The forces and 
responses are related by a mathematical operator. 
Block programming starts with the selection of basic 
definition and rules. 


MATHEMATICAL OPERATOR (0) 


The operator, or combination of operators, describes 
the functional relationship between quantities. The 
most common operators are the following: 

Summer 


Constant Multiplier (potentiometer, commonly 
referred to as "pot") 


Inverter - -1 Sign Changer 
1 
Integrator - > (La Place Transform Notation): 
Differentiator - s (This operator is never used in 
an analog simulation if itcan 
be avoided, as it usually can). 
Variable Multiplier 
Arbitrary Function - f(x) 


Analytical Function Trigonometric, Log, 


Hyperbolic, etc. 


Inverse Operator - (0)? S, > are inverses 


ELEMENT 


The element is the basic unit of the block representation 
of asystem. It hasa pair ofterminals or nodes, with 
a value associated with eachnode. It alsohas a trans- 
ference quantity between nodes, or ''through" the ele- 
ment which is functionally related by some operator to 
the node-pair value. The transference quantity may 
also be referred to as the branch transference or trans- 
mission. 


The elements required to depict the passive lumped- 
parameters of most physicalsystems are very few in 
number. Inparticular, five basic elements will suffice 
for many linear electrical, mechanical, and thermal 
systems. These aretheSummer, Constant Multiplier, 
Sign Changer, Integrator, and Differentiator. 


BLOCK 


A block is the computer mechanization related toa 
physicalelement. The input/output variables of acom- 
puter block are voltages andcurrents, but the operator 
relationship is the same as the simulated element. 


Table I gives examples of single elements and their 
corresponding blocks. Table IIgives afew commonly 
occurring two-element combinations. Combinations 
of more than two elements can be generated as shown 
in the following examples. Usefultables of more com- 
plicated blocks or transfer-function simulations appear 
in many publications. 


DRIVING SOURCE 


A driving source is a source of energy or power for a 
network of elements. 


SYSTEM STABILITY 


A stable system is one whose responses are bounded 
(finite amplitude limit) for any finite input. An un- 
stable system is one whose responses are not bounded. 


NETWORK 


A network is an interconnected set of elements and 
sources. Elements may be joined together at their 
nodes, providing the connected nodes have common 
dimensions and value, and provided the transference 
quantities have common dimensions. 


This discussion concerns elements and networks in 
which the following two rules apply: (Kirchoff's laws) 


1. Thealgebraic sum of the node values around 
any closed path (loop) in a network is zero. 


2. Thealgebraic sum of the transmission quan- 
tities at any node is zero. 


: Korn and Huskey, "Computer Handbook", Chapter 2, 
McGraw-Hill, New York, 1962. 


TABLE I. Single-element Computing Blocks 


PHYSICAL SCHEMATIC MATHEMATICAL 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 
Summer > X | _ | >— 
Y=- (X, + Xo) Xo 


Constant 
Multiplier 


(POT) 
Y = KX 


Differentiator 
Y =sxX 


Linear 
Spring 


Torsional 
Spring 


Heat 
Conductor 


Viscous 
Damper 


TABLE I (Continued) 


PHYSICAL SCHEMATIC MATHEMATICAL 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 


Electrical 
Capacitor 


Viscous 
Damper 


Mechanical 
Mass 


Thermal 
Mass 


Electrical 
Inductor 


Mechanical 
Mass 


COMPUTER PROGRAM MECHANIZATION a) assigning aposition inthe mechanization 
COMMENTS diagram to all necessary variables and 
operators. 


b) accomplishing the necessary intercon- 


1. Operational amplifiers operate with the junc- nections so that each output has its pro- 
tion at virtually zero potential; the currents per inputs. 
into the junction through all the input and 
feedback paths must sum to zero; and there c) checking to see that all input/output rela- 
is asign inversion between the amplifier in- tionships in the mechanization are ‘satis- 
put and output quantities. fied. 
2. Anyclosed-loop ina problem mechanization EXAMPLES OF NETWORK PROGRAMMING 
simulating astable physical system will gen- 
erally have an odd number of sign inversions. The following examples illustrate some of the tech- 
niques used in applying the block programming meth- 
3. In mechanizing an equation, set the highest ods. Examples are given for linear elements only. 
order derivative only on the left of the equals Non-linearities can be included simply by replacing 
sign to obviate the necessity for differen- the potentiometer corresponding tothe non-linear ele- 
tiation. ment with a multiplier, function generator, or other 
suitable computing element. Backlash, hysteresis, 
4. The mechanization is accomplished by: stiction, etc can be included as required. 


TABLE Il. Two-Element Computing Blocks 


MATHEMATICAL 
SCHEMATIC RELATIONSHIP PROGRAM BLOCK 


MASS - SPRING - DAMPER SYSTEM 


Xo1Vo 


MOTION IN X 


F DIRECTION ONLY 


Element - force equations 


= K 


Spring K, : F. 1 


Spring K,: F_, = K Initial values att = 0. 


52 XK, - X 


g | 2 1) 
X (0), X, (0), vy (0), vy (0) 


Damper B: F,. = B (v, - v,) 


D 1 


Forces acting on masses 


F = -Foy t Feo D 


Fu2 = “Fs 7 Fp + F 


Steps to mechanize 


1. Assume Xi -Xp, “V4 Vo and F are available 


2. Draw blocks for M,> My 


3. Draw blocks for K,; Ky, B 


4. Make input connections to mass blocks to produce xX, “Xo, -V 


——--—— —  — 7 - 
Pema | tt pe 


| 
| 
| 
| 
| 
|" 
| 
| 
| 
| 
| 


SPRING! | MASS! !|DAMPER | SPRING2! MASS 2 


——-L__ tee LLL LL 


The above mechanization makes each parameter avail- 
able as an individual adjustment. The programming 
is easily accomplished by generating the forces im- 
posed on the elements and then operating on the forces 
to obtain velocity and displacement. Any or all of the 
initial conditions, K,(0), X2(0), vy(0) and v9(0), may 
be zero, or may also be variable. 


Note that the block representing the damper contains 


two summers. The first summer produces (vy - Vg), 
and the Second summer acts only as a signchanger. 
The two-amplifier, one-pot combination can be reduced 
to a more simple two-pot combination by connecting 
V1 and vo through B pots to the mass block inputs. 
The less complex alternate isshown below. This con- 
figuration requires fewer amplifiers, but variations 
in M,; Mo, or B necessitate two or three adjustments. 


* Do not connect to “Vy to get pot setting of B instead 


_.° A coefficient pot output cannot be connected 
1 


to another coefficient pot input. This is a practical 
limitation of the computer hardware. 


Rigid Bar Supported on Springs 


e 22 1 1 F 
Q = 
sI 
FH = -K, Y 
Fo = -Ky Yo 


Rigid bar of mass M, moment of inertia about C.G. of 
I. Assume 6smallsuchthat 6+singo. Vertical motion 
only. Yo is displacement of C.G. 


: - FL + Fa + F 
0-7 sM 

Yi = Yo 4,9 

Yo = Yo + £99 
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One-Dimensional Heat Transfer 
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M-MASS 


C-SPECIFIC HEAT 


HERMAL 
MASS 
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AL 
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Electrical R-L-C Networks 


RESISTOR (CAPACITOR) RESISTOR CAPACITOR j NOUCTOR RESISTOR 


Lot dL LL 


Io Iz 
1, I, 
Fi - Be = sc, Ba B3 > 8c, Ez = 1, R, 
E 
I = — lL =L-!I 


TWIN-T NETWORK 


; _ Fy” Es - _ Igo te 
2 R, 3 SC, 

Pe ier: 
3 R, 

-(I, +I,) 

g tly. 

(E, - Ey) = sC, 

I, - 1, 


E, = (I, - 1, +1, +1,) Ry 


E, = (Ey - Ey) +E4; E, = Ey - (Ey - Ep) 


-(T,-Ia+Iz+Ig) 


CRs 
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SCALING 


The mechanization diagrams associated with the pre- 
ceding examples show the computing element inter- 
connections without indicating potentiometer settings, 
amplifier input resistor values, andintegrator capac- 
itor values. These are unscaled mechanizations. 
Scaling is required to relate the computer voltages to 
the physical problem variables. 


Tne analog computer is voltage limited to a practical 
range of operation of +100 volts as a maximum and +1 
volt or zero aS a minimum. It is usually not possible 
tonumerically equate analog voltages to physical vari- 
ables. That is, the magnitude of physical variables 
must be scaled to fall within the useful range of the 
machine. This procedure is called amplitude scaling. 


The analog computer is also limited as tothe speed with 
Which it will solve a problem. Practical bounds on 
solution time are 100 seconds as a maximum and 100 
milliseconds as a minimum. Events in the physical 
world usually occur in time intervals which fall out- 
side these limits. Thus, the simulation of real world 
phenomena with the analog computer is ordinarily 
faster or slower. Theprocedure for relating computer 
time to physical time is called time scaling. 


Correct scaling is animportant factor in reducing sim- 
ulation errors. An error analysis of any particular 
program is about as complex as the problem which is 
being solved by the program. Consequently, except 
for unusual cases, an analysis of error is not made. 
An estimate of the accuracy of the program can be 
found from check cases for whichthe answer is already 
known. Error is reduced by scaling the program so 
that all potentiometer values and amplifier gains are 
reasonable. Frequently, in the process of scaling a 
problem, it will be found that some parts of the model 
are not significant and can be eliminated. Scaling is 
a good check on reasonability. 


Scaling is anart, not a science. A good deal of exper- 
ience is required to become proficient. Scaling canbe 
done inseveral ways. Thetechniques discussed in this 
chapter are those in most prevalent use among analog 
programmers. 


If X is a physical variable corresponding to computer 
voltage, V, then a scale factor, a, is chosen so that 


Ve=a-xX 


will fall within the practical range of the computer. It 
is good practice to make V as large as possible. If 
is physical (real) time corresponding tot,, the com- 
puting time, then these two can be related by a scale 
factor, N: 
m 


t = Nt, at” = N™at™ 
Cc p’ c p 


It follows that 
fcFat, = Nf (Foat,. (1) 


N is chosen so that the problem solution time on the 
analog computer falls within reasonable limits. Prac- 
tically, this is accomplished by a choice of N which 
will not result in prohibitive amplifier gains. 


The general procedure for scaling is: 


1. Generate anunscaled program having a math- 
ematical structure which agrees with the pro- 
blem (model) to be solved. 


2. Identify the location of all voltages corres- 
ponding to physical variables. 


3. Associate a scale factor with each of these 
voltages (e.g. Vy = aX). 


4, Label the program with the scaled physical 
variables (i.e. aX) rather than the voltages. 


5. Estimate maximum values of the physical 
variables. 


6. Choose the amplitude scale factors andthe 
time scale factor, N, so that all pot settings 
and gains are reasonable. 

In order to accomplish step (4) it is necessary to know 
how input and output scale factors are related for 
various computer elements. 
POTENTIOMETER 

Physical equation: Y = bX 


Scale factors 


y ¥y 
4y 
Scaled equation vy = Pa V,. 
Program aX ( ) a,¥ 
ba 
a 
x 


Note that a scale change can be made withapot. Assum- 
ing either a, or has been established previously, 
the other is chosen so that ba,/ay is a reasonable pot 
setting. The pot can be used solely for a scale change. 
In that case b = 1. 


SUMMER 


Physical 
equation: Z = X + Y 


Scale 
factors V. =a Z, V_=a X, V_= ayy 


Scaled a, a, 
equation: V_ =— V, + — V 

Z a a 

x y 
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Now, if only one-gain inputs for the summer are used, 
then it must be true that a =a =a and the pro- 
gram is y 


-a,X 


Q7Z 
—ayY 


Thus, it is seen that input scale factors must be the 
same for a one-gain summer and that a change of 
scale cannot be made through asummer. However, 
a different input scale factor can be used fora 10-gain 
input. The program for this case is, 


-a,X 


_ a,Z 
ayY lO 


where it must be true that a. = 10a, . 


INTEGRATOR 


As with the summer, the input scale factors (except 
for the 10-gain input) must be the same. However, a 
change of scale can be made with an integrator. Both 
amplitude and time scaling can be done. Scaling for 
an integrator is derived below, where C is the value 
of the capacitor in yufd, and R is in megohms. 


Physical 
equation: Y =f Xdt, 


Scale 
factors Vy= a,¥; Vy. = aX, th = Nt, 
Scaled a, 1. 
equation: V, = a ON fv, at, 
Thus, the integrator gain, aa must be 
a. fy 
RC ay (2) 
and the program is 
aX NayY 
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MULTIPLIER 
Physical 
equation: XY = 2 
Scale 
factors: V. = aX, V_=a_ Y, V_=aZ 
x x y y Z Zz 
Scaled 
equation: V_V 100a 
xX y Z _y 
100 aa Z 
x y 
Program: 
OxX a,Z 
ayy 
It must be true that 
100a 
4 =1, 
aa 
x y 


From this requirement itis seen that the relationship 
among the scale factors is the same as thatfor input and 
output multiplier voltages: 


aa 
a, = y 
100 
DIVIDER 
Physical 


equation: = = Z 


Scale 
factors 


Scaled V 


‘on: xX _4Y 2. 
equation: OF . 100a,, V, 


Scale 
factor a. 
constraint: a = 100— 
Z a 
y 
Program: 
aX a,Z 


SQUARE ROOT 7 
OO Estimates: (ee =8x10 3. Xmax. = 1 
max 


Physical p 


equation: Y = /X 


The unscaled program is: 


Scale 
fact = _ 
ors vy a¥, V,. ax 
Scaled a 
equation: Vy > y . 10 /V, 
10 /a., 
Scale 
factor 


constraint: = 
raint a, 10 fa, 


Program: 


The scaled program is obtained by assigning literal 
scale factors, determining their numerical values, and 


FUNCTION GENERATOR calculating the coefficient potentiometer settings. 


Physical 
equation: Y = f(x) 


TTT -100 
Scale 
facto = = 
rs Vy ats Vy aX 

Scaled 
equation: V_= af(— 

y y a, 

EXAMPLE: 


Mass-spring-damper system 


cKX cK 


N@M N2 aM 
1 1 
The = and 57 factors appear 


due to the inherent gain of N 
through each integrator, from 


Eq. (1). 


In general, the determination of the numerical values 


ayer Spring F F = KX involves atrial and error approach. Any convenient 
q * Spring soree, ts parameter can be picked for the starting point. Ifa 
ax selected value leads to unreasonable potentiometer 
Damper Force, Fy = Ba settings or unreasonable amplifier gains, then new 
p values may be required. Reasonable pot gains are 
between .1 and1. Reasonable amplifier gains are 
Mass q2 -(F_ + F) between 1 and 20. Amplifier gains should be integer 
Acceleration, —3 = —a values. Non-integer values can be obtained from a 
dat potentiometer, in combination with gain-of-10 inputs, 
if necessary. 
Be -4. -5 

M=1;B=4x10 ;K=6.4x 10 1. Determination of a. In order to utilize the 
ax full dynamic range of integrator 2, set 

at (0) =0, X (0) = .95 aX max = 100 (volts). 

p 
a = 100 
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Potentiometer 1 setting. 


= .95 (3) 


Determination of a At this time, a trial 


value for N maybe selected. Analog solution 
times of .1 to 100 seconds, or computer nat- 
ural frequencies from .1 to 100 radians per 
second are the preferable operating ranges. 
Anestimate of the physical time or frequency 
ranges can be used to Select a value for N. 


For this example, a valueof N = .01 appears 
reasonable. 


b — 100 (volts) 


N (dX/dt max. (4) 


Looking ahead, note that integrator 2 will 
aN 
5 
herent in the integrator from Eq. 1. There- 


require a gain of The gain of N is in- 


fore, the value for the integrator will be 


= It is highly desirable that this factor be 


aninteger, preferably either 1 or 10. There- 
fore, b = 100 will be used. The effect of this 
is to force integrator 1 to utilize something 
less than its full dynamic range, but the effect 
is not serious inthis example. As a general 
rule, it is not possible to scale so that all 
amplifiers work over their full range. 


4. 


Potentiometer 2 setting; 


P2 = x - (6.4x 107° 
N aM 


)e (5) 


The value of c is chosento assure that a rea- 
sonable value for P2 is obtained and to assure 


that — is an integer. 
c = 100 will satisfy both conditions. 


P2 = .64 


Potentiometer 3 setting: 


_ cB _ 
P3 = yom . 04 


This setting for potentiometer 3 is below the 
desired minimum potentiometer setting value 
of .1. Inthis example, however, it is not pos- 
sible to have a large setting for potentiometer 3. 
This follows from a consideration of the phy- 
sical problem. The damper, B, was chosen 
tohavea very slight effect on the mass-spring 
system, so it necessarily follows that the 
portion of the analog program related to the 
damper willhave onlya slight effect. It should 
be noted that abnormally low potentiometer 
settings will frequently arise in simulations of 
systems containing negligibly small elements. 
It should be recognized that the low settings 
stem from the actual system characteristics. 
The programmer, therefore, need not spend 
time in futile attempts to improve the scaling. 


The final program with numerical values is: 


dX 
——=|0"4 
dtp Mt 


2 
d“xX - 


p 


GENERAL COMMENTS ON TIME SCALING: 


Tne following points concerning time scaling are of 
varticular importance: 


1. 


If all integrator gains are simultaneously 
changed bya factor k, the time scale changes 
from N tor, and all computer frequencies 
increase by k. The amplitude scaling, how- 
ever, does not change. All numerical scale 
factors, all potentiometer settings, and all 
computer problem voltages remain the same. 
This property of amplitude invariance under 
time scale change allows the programmer to 
easily speed up or slow down the computer 
solution. 


The most convenient way of changing all inte- 
grator gains simultaneously is by changing 
the capacitors by a fixed factor. Integrator 
capacitors can be simultaneously changed by 
means of the computer Time Scale control+. 


An event which occurs in a fixed real-time 
interval can occur in several different com- 
puter time intervals depending on the choice 
of integrator capacitors, that is, depending 
on the time scale factor. 


Events of identical character which occur in 
different real-time intervals can all be made 
to occur in the same computer time interval 
by appropriate time scaling. 
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on the SD 10/20 and SD 40/80 analog computers, a master Time Scale switch is a standard feature. This per- 
mits fast and convenient switching of capacitors that are patched -- up to 1000:1 time scale (x10, x100, x 1000, 
depending on patch panel connections). 
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